Microscopic correlation functions for the QCD Dirac operator with chemical potential 
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A chiral random matrix model with complex eigenvalues is solved as an effective model for QCD 
with non-vanishing chemical potential. We derive new matrix model correlation functions which 
predict the local fluctuations of complex Dirac operator eigenvalues at zero virtuality. The parameter 
which measures the non-Hermiticity of the Dirac matrix is identified with the chemical potential. 
In the phase with broken chiral symmetry all spectral correlations of the Dirac eigenvalues are 
calculated as functions of quark flavors and chemical potential. The derivation uses the orthogonality 
of the Laguerre polynomials in the complex plane. Explicit results are given for any finite matrix 
size N as well in the large- iV limit for weak and strong non-Hermiticity. 
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Random Matrix Models (RMM) have been a useful 
tool in Theoretical Physics for a long time. In many 
physical systems the local fluctuation properties for ex- 
ample of the energy levels are universal and can be suc- 
cessfully described by RMM, where we refer to § for a 
review. Although in the generic situation the Hamilto- 
nian and other physical observables arc Hermitian thus 
having real eigenvalues, there exist also important cases 
where complex eigenvalues occur. As examples we men- 
tion localization in superconductors dissipation and 
scattering in Quantum Chaos || or Quantum Chromo- 
dynamics (QCD) with chemical potential 

Much less is known so far for spectral correlations of 
complex eigenvalues derived from RMM. Although the 
first results date back to Ginibre || where the correla- 
tions for the complex Unitary Ensemble labeled by the 
Dyson index /? = 2 were calculated, progress has been 
slow. The correlation functions of the ensemble with real 
non-symmetric matrices (J3 = 1) are still unknown. Re- 
sults for quaternion matrices (fl = 4) were obtained more 
recently in || Q and the inclusion of Dirac mass terms for 
(3 = 2 in ]8|. Furthermore, it has been realized by works 
of Fyodorov and collaborators || that different regimes 
of complex eigenvalues exist, the weak and strong non- 
Hermiticity limit. In the present work we wish to extend 
the knowledge about complex matrix models to the so- 
called chiral models. 

Chiral RMM of real eigenvalues have been introduced 
to describe the local fluctuation properties of the Dirac 
operator in QCD at the origin jlO| . The low energy spec- 
trum of the QCD Dirac operator is a very sensitive tool to 
study the phenomenon of chiral symmetry breaking [ ll| . 
The predictions of the different chiral RMM ensembles 
have been very successful in describing the dependence 
on the gauge group and its representation, the number 
of quark flavors and masses and the topology j^j. In 
particular the topology dependence has been very useful 
in comparison with new developments in lattice gauge 
theory, admitting to incorporate an exact chiral symme- 
try p3[ . By now the field theoretic origin of the RMM 
description has also been well understood H]. 

On the other hand lattice simulations in the presence 
of a chemical potential ii, which renders the Dirac eigen- 



values complex, remain extremely difficult, as reviewed 
in |^| . Although recent progress has been made 
along the phase transition line the general phase diagram 
remains unexplored for fi 0. It is in this context that 
analytical knowledge of microscopic correlation functions 
from complex chiral RMM could be very useful, in view 
of its success in predicting real eigenvalue correlations. 

Chiral RMM including a chemical potential have been 
already studied in several works J§ [H| In || 
the nature of the quenched limit has been analyzed in 
such a model, and the global density of complex eigen- 
values together with its boundary have been calculated 
as a function of /i. The phase diagram of QCD in the 
temperature density plane has also been predicted with 
such a model JlS| . Recently complex Dirac eigenvalues 
calculated on the lattice have been confronted to a com- 
plex matrix model on the microscopic scale given by the 
inverse volume in the bulk of the spectrum JnJ. The 
nearest neighbor distribution along a given direction in 
the complex plane was considered and a transition from 
the Unitary to the Ginibre ensemble was observed at in- 
creasing fi, ending in a Poisson distribution. 

Our aim here is to provide more detailed information 
by calculating all microscopic correlation functions in the 
complex plane, both in the limit of weak and strong non- 
Hermiticity. Our model gives a natural chiral extension of 
the ensemble treated in [|[ |) . Although our main moti- 
vation is the application to QCD lattice calculations with 
/i =/= the new correlation functions we derive may find 
other applications as well, such as in the fractional Quan- 
tum Hall effect pOl or two-dimensional charged plasmas 
pit . We restrict ourselves to work with a Gaussian chiral 
RMM. In particular we will not touch the issue of uni- 
versality here although we expect that in analogy to Q 
the same correlation functions hold for a more general 
weight function at large- A?\ 

We start our investigation by defining our model. The 
chiral random matrix partition function in terms of the 
complex eigenvalues Zj=i,...,jv of a complex N x N matrix 
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T-r-A is defined as 
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where we have introduced the weight function 
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It contains Nf massless quark flavors in a fixed sector 
of topological charge v, where a = Nf + v. We have 
taken the absolute value of the Dirac determinant but 
the weight depends on the real and imaginary part of 
the eigenvalues. The weight has been chosen to be Gaus- 
sian with variance (1 + t)/2N for both the Hermitian 
and anti-Hcrmitian part of J, H and A respectively. 
The parameter r € [0, 1] controls the degree of non- 
Hermiticity and will be related to the chemical poten- 
tial \x in QCD below. The Jacobi determinant from the 
diagonalization of J j^] yields the Vandermonde determi- 
nant Ajv(zi, zn) — n^>i( z fc~ 2; 0- The reason for the 
model eq. (0) being chiral can be seen as follows. For 
real eigenvalues the chiral ensemble [[To] is usually de- 
fined on the real positive line, J °° <iAA a exp (— NX). By 
substituting A = y 2 it can be mapped to the full real axis, 
jToc dy |y| 2o+1 exp (—Ay 2 ), where y corresponds to a real 
eigenvalue of the Dirac operator p |2j| . In the latter pic- 
ture the continuation into the complex plane is straight- 
forward, leading to the ensemble (|l ) . The only difference 
to the massless non-chiral ensemble || is the additional 
power in \z\. Furthermore, the orthogonal polynomials 
for the ensemble (|l|) are given by Laguerre polynomials 
L^(z 2 ) in the complex plane. 

A chemical potential fi is included in the QCD action 
by adding /X70 to the Dirac operator p and it renders 
its eigenvalues complex. Our aim is to describe the lo- 
cal fluctuation of small eigenvalues close to the origin, 
because of their importance for chiral symmetry break- 
ing through the Banks-Casher relation p4| . For fi not to 
completely dominate the Dirac determinant we will also 
restrict ourselves to small values of \i. 

In Q a slightly different matrix model was introduced 
replacing p by a chiral random matrix and keeping the 
additional term /170 explicitly. Here, we treat p + /i7o 
as a chiral, complex matrix instead. The advantage of 
our model is that we can directly calculate all micro- 
scopic correlation functions. We will take the model 
Q to relate our parameter r to \i by comparing to the 
macroscopic spectral density p(z) and its boundary cal- 
culated there as a function of fi. In the limit of small 
/i the spectral density of Q becomes approximately con- 
stant, p(z) = l/47r/x 2 , and it is bounded by an ellipse, 
x 2 /4 + y 2 /An 2 = 1, where z = x + iy. This behavior 
can also be observed for lattice data with small fj, (e.g. 
in ]l9| , p5|). The macroscopic density in our model can 
be read off from |E6| since the Dirac determinants are 



subdominant in the macroscopic large- N limit: 
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We therefore identify 

V = (1-r 2 ) 



< 1 
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valid for small chemical potential and r close to unity 
meaning small non-Hermiticity. For large values of // the 
eigenvalue density on the lattice is no longer constant 
and develops a hole in the middle (see e.g. 1 19 , ^5|). We 
will also see such a hole develop in the microscopic cor- 
relations in the limit of strong non-Hermiticity, as shown 
in Fig. § below. 

After having identified all parameters in our model eq. 
(jl]) we turn to its solution using the powerful method of 
orthogonal polynomials || . We only state the results and 
refer to ]27j for details. All eigenvalue correlation func- 
tions are first given for a finite number of eigenvalues 
N and then in two different large- TV limits correspond- 
ing to weak and strong non-Hermiticity. The orthogonal 
polynomials in the complex plane are defined as 



dzdz* w {a) {z)P y k a >{z)Pl a >{z*) = 5, 
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Following standard techniques || the knowledge of the 
kernel of orthogonal polynomials 



N-x 



(6) 



fe=0 



allows to calculate all fc-point correlation functions 

■ (7) 
is given 
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The result for the orthogonal polynomials eq. 
in terms of Laguerre polynomials 
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with the normalization integral 
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and P a+ |(a;) being the Legendre function. All fc-point 
correlation functions then follow by inserting the polyno- 
mials into eqs. (^) and (0). In our results the parameter 
a = N f + v can be kept real (with a > —1). For example 
we can set a = — 5 as a check, recovering the even subset 
of the Hermite polynomials in the complex plane [E0[. 
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Since we did not find eq. (g) in the literature we briefly 
sketch its derivation. Performing a change of variables 
z — > e %<p z in the normalization integral eq. ([)]) we obtain 

i /Nz 2 ' 



1 = 



with 



exp 
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r 2 (l-e 2 ^) 

= (1 _ T 2 e 2^) 



(10) 



(11) 



and the average taken with respect to w^ a '{z). Multi- 
plying both sides of eq. (|l^) with ((1 - u)(l - u*)) _a_1 
and recognizing the generating functional of the Laguerre 
polynomials we obtain the desired orthogonality relation, 
given properly normalized in eq. (^). 

After giving the exact solution for finite N we turn 
to the large- N limit. We first consider the weak non- 
Hermiticity limit. Following Q we take the limit t — > 1 
such that the combination 



lim N(1-t 2 ) = a 1 = ANi? 



N 



(12) 



is kept fixed. Because of the identification eq. (Q) we 
consequently also rescale /i going to zero when N — > oo. 
In other words the weak non-Hermiticity parameter a 2 
directly measures the chemical potential in the micro- 
scopic scaling limit. Such a rescaling is similar to that of 
the quark masses . It has been already mentioned in 
|L7| that in a RMM the numerical effort to obtain con- 
vergence grows exponentially with 7V/i 2 . Keeping it fixed 
here should make a comparison to data feasible. 

Furthermore, we also rescale the complex eigenvalues 
keeping 



iV(5ie z + i3?m z) = Nz = £ 



(13) 



fixed. The matrix size N corresponds to the volume on 
the lattice. This defines our microscopic origin scaling 
limit in the complex plane. The kernel eq. (o) and cor- 
relators eq. (fjj) also have to be rescaled with the mean 
level spacing 1/N of the eigenvalues. 

In order to obtain the microscopic kernel from eq. @ 
we replace the sum by an integral, J2k=o ~* ^ Jo 
where t — 4r, and use the asymptotic limit of the La- 
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guerre polynomials to finally arrive at 
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J a (V2t^)J a (V2tC 2 ). (14) 



The microscopic, weakly non-Hcrmitian correlation func- 
tions obtained from eq. (Q) are our first main result: 
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As an important check in the Hermitian limit a 2 
corresponding to t = 1 the universal correlations jlOj |22|] 
of the chiral RMM with real eigenvalues are reproduced. 

To give an example for eq. (^) we have depicted 
the quenched microscopic density in part of the complex 
plane in Fig [j]. The other directions follow from sym- 
metry. The oscillations known from the real case |lC| ] 




FIG. 1: The quenched microscopic density for a 



0.6 



spread into the complex plane, indicating the locations 
of the individual eigenvalues. 

We now turn to the strong non-Hermiticity limit. In 
this limit r 6 [0, 1) and consequently also /i from eq. (^) 
is kept fixed in the large- TV limit. The eigenvalues are 
now rescaled with the square root of the volume ^ [| , 



y/W (3?e z + i$tm z) = Vn z = £ 



(16) 



defining our microscopic origin limit at strong non- 
Hermiticity. In this limit the infinite sum in the kernel 
eq. (^|) can be evaluated using standard formulas for La- 
guerre polynomials and we obtain 
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The correlation functions in the strong limit then read 
k 
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det 
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(18) 

where we have suppressed the normalization constant. 
They can also be obtained from the correlators eq. (|l5| ) 
in the weak limit by taking the limit a — > oo there and 
identifying a 2 = 1 — t 2 . They differ from ref. 0] due to 
the different interaction term in eq. ([!]). As an example 
for eq. fllq ) we give the quenched microscopic density 
in Fig. H The microscopic spectral density develops a 
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FIG. 2: The quenched microscopic density for r = 0.5 

hole at the origin and becomes flat for large values. This 
has indeed been observed in quenched lattice data for in- 



termediate values of the chemical potential ^5|. At 
large fj, the eigenvalues form a thin ring. The qualitative 
picture of the microscopic density in Fig. ^ remains un- 
changed in the presence of massless flavors Nf, although 
the barrier present in the figure gets flattened out. 

In summary we have introduced a new class of a chiral 
RMM having complex eigenvalues. All fc-point correla- 
tion functions have been calculated explicitly at finite- 
TV as well as in the limits of weak and strong non- 
Hermiticity. The parameter r that governs the non- 
Hermiticity has been related to the chemical potential 
fi as it occurs in the QCD Dirac operator. While the 
microscopic density shows qualitative features of Dirac 
operator eigenvalues calculated in lattice QCD a quanti- 
tative comparison with data remains to be done. 
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